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1 Introduction 

The analysis of the polarization characteristics displayed by optical systems can 
be- divided' into two categories: geometrical and physical. Geometrical analysis 
calculates the change in polarization of a wavefront between P U P“ S ‘ n an opti- 
cal instrument [1,2, 3, 4]. Physical analysis propagates the polarized fields wherever 
the geometrical analysis is not valid, i.e. near the edges of stops near images, in 
anisotropic media, etc. [1,5, 6, 7, 8]. The changes, geometrical and physical, po anza- 
tion causes in the performance of lens and mirror systems are readily calculated by 
several commercial computer codes [9,10,11,12], The inverse problem of designing 
a system with specified polarization characteristics is more difficult. Examples l 
the literature include a polarization compensated polarizing microscope[13] and a 
telescope with ultra-low polarization for a solar polarimeter [14]. Design requires a 
fundamental understanding of the origin of polarization aberrations an ow ey 
change with both the optical and coating prescriptions for a system. Polarization 
aberration theory provides a starting point for geometrical design and facilitates 

subsequent optimization. . ., . .. _ 

Chipman has derived several polarization aberration expansions similar to the 

classical wavefront expansion for rotationally symmetric systems valid for weak : sec- 
ond order aberrations [15,16,2], Reference [17] explores polarization aberrations 
graphically using a “symmetrized” second order expansion valid for strong aberra- 
tions but does not describe a method to calculate the aberration coefficients, n 
this paper, we calculate and discuss an exponential expansion of the polarization 
aberrations valid for strong polarization aberrations through fourth order . The 
results are applied to the interpretation of polarization raytracing results. 

The polarization aberrations described in this paper arise from differences m the 
transmitted (or reflected) amplitudes and phases at interfaces. In contrast, classical 
wavefront aberrations arise from differences in optical path length as rays propagate 
between interfaces [18]. Figure 1 shows the calculation ot the optical path length 
W and the optical path length and the polarization J along a ray. Repetition o 
the calculation depicted in Figure 1 (a) for multiple rays and wavelengths samples 
the wavefront aberration function. Repetition of the calculation 1 (b) for multiple 
rays and wavelengths samples the wavefront aberration function and the polariza- 
tion aberration matrix (PAM) of the system. The PAM describes the variation in 
polarization with object coordin. pupil coordinates, and wavelength. This P a P er 
calculates the PAM for isotropic rotationally symmetric systems through fourth 
order and includes the interface phase, amplitude, linear diattenuation (defined 
in Table 1), and linear retardance aberrations. Polarization aberrations resulting 
from propagation through anisotropic media such as crystals are not considered in 
this paper. For propagation through anisotropic crystals, the propagation terms 

expfi'Wi in Figure 1 (b) would be replaced by Jones matrices J nm . 

The order of an aberration term referred to in this paper is the order of the 
wavefront representation (n), not the order of the transverse aberrations (n - 1). 
Thus, defocus and tilt are second order aberrations, while spherical aberration, 
coma, and astigmatism are fourth order aberrations (not third order aberrations). 

Section 2 discusses the exponential form of Jones matrices used in this paper. 
Section 3 introduces the PAM in Jones matrix form. In Section 4, the exact cal- 
culation of polarization aberrations through polarization ray tracing is described. 
Section 5 presents the coordinate system used in this paper. Section 6 discusses the 
paraxial approximation including: thd pafhjcial PAM for f single surface, paraxi 
angle of incidence, and the paraxial orientation of the plane of incidence. In Sec- 
tion 7, a Taylor series simplifies coating dependence of the single surface paraxial 

lThe preliminary version of this paper used did not use the exponential form which resulted in 
increased complexity of computation and interpretation. 
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Table 1: Polarization terminology 


*[trm Definition 

Diattenuation The prope7ty“^f having an intensity transmittance 

which is dependent on the incident polarization 

state 

Polarization The property of altering the polarization state of 

light. Polarization includes the subsets of diattenu- 
ation and retardance. 

Polarizer An optical element which transmits a fixed polariza- 

tion state independent of the incident polarization 
state. Examples include dichroic sheets (Polaroid) 
and Glan-Thompson prisms. 

Polarization element Any optical element showing polarization. Exam- 
ples include retarders, polarizers, and metallic in- 

t€ r f* aces 

Retardance The property of having a phase or optical path 

length which is dependent on the incident polariza- 
tion state^ ___ - — 

Note: There is a distinction between polarizer and polarization element. Chipman 
provides a more detailed discussion [2]. 


PAM Section 8 calculates the paraxial PAM for a system of isotropic rotationa y 
symmetric elements through fourth order. A general discussion of the terms is 
contained in Section 9. Section 10 contains a detailed discussion of the vector a er- 
rations (defined later) comparing and contrasting them with classical scalar p ase 
aberrations. Section 11 discusses interpretation of polarization raytracing in the 
• context of the aberration theory results. Appendix A contains paraxial expressions 
for polarization basis vectors. Appendix B examines the polarization by uncoated 
interfaces. Appendix C lists the polarization aberration coefficients for a system ot 
isotropic rotationally symmetric elements. 


2 Jones Matrices 

In this section, we present the Jones matrix formalism[19,20] for the analysis of 
polarization as used in^this paper. 

The Jones vector U{t) is 



where U q {t) and U r {t) ar» the projections of the electromagnetic field on any two 
orthogonal basis states q and r. The Jones matrix J relates the incident U and 
transmitted fields U f 

O' = iU ( 2 ) 


and is a 2x2 matrix with complex elements 
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Table 3: Physical significance of the exponential polarization coefficients 


Coefficient 

Matrix 

Physical significance 

ClQ 

<7q 

3t(a„): 

Polarization independent amplitude 



3(a 0 ): 

Polarization independent phase 

3i 


*(oi): 

Linear diattenuation along the coordinate 




axes 



0(ai): 

Linear retardance along the coordinate 




axes 

02 

<*■2 

3J(a 2 ): 

Linear diattenuation at 45 degrees to the 




coordinate axes 



0(a 2 ): 

Linear retardance at 45 degrees to the co- 




ordinate axes 

33 

<T3 

3 ?(a 3 >: 

Circular diattenuation 



9(a 3 ): 

Circular retardance 


U i p (h, p, A) = Jjyi(h, P< A )U n p(h , p, A) 
where the Jones matrix is 


(5) 


^ syt(h , pt A) 


f Jjy»,n(h| p, A) jjyj,12(h,P, A) \ 
y i»y*, 2 l(h, Pi A) Jjyj,22(h, p, A) J 
expfajy^otf'O + ajyj.lO’l + a jyj : cr 2 + a aya 


( 6 ) 


and h is the object coordinate, p is the pupil coordinate, and A is the wavelength. 

It is convenient to separate the Jones matrix for an imaging system into a polar- 
ization aberration matrix (PAM) which describes the aberrations and a "quadratic 
phase” characteristic of ideal imaging systems 


J JVJ (h,p,A) = J(/», P, A) exp 



(7) 


where 3(h,p, A) is the PAM, / is focal length, and k = 2a-/A is the wavenumber. 
Both the wavefront aberrations and the finite extent of lens are described by the 
elements of J (the finite size of the lens is an “aberration” which reduces resolution 
from that predicted by geometrical optics.) In the limit of a non-polartzing optical 
system, the PAM has the form 

J(A,p, A) = P(h, p) exp [jkW{h,p, A)] q l ) ^ 


where P{h,p) is the pupil function which describes the amplitude transmittance 
of the pupil and W(h,p,X) is the wavefront aberration function. In succeeding 
sections, quasi-monochromatic light is assumed and the explicit wavelength depen- 
dence A is suppressed. 


4 Exact Polarization Raytracing 

In this section, we describe the procedure for calculating the change in polariza- 
tion along a ray through a system of isotropic surfaces. " The technique is called 
polarization raytracing [ 23 , 24 ]. ... , 9 ^ ; i' V*V. . _ 
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J ,(i,A) = . 

( <*»»(**) c08 a »«(*«) sin ^ ) . ( 19 ) 

~ \ -a p ,(»,) sin 9 q a P1 (i q ) cos 9 <, ) 

Equation (19) is exact Jones matrix for an isotropic interface- 

do not display circular polarization or circular retardance. I a * 1 h ff 

on a system of surfaces, the polarization along ft is found by cascading the effects 

from each surface 


J(/?l) — 3q(iq,9 q ) 


m 


The set of all Jones matrices for a system for each possible ray P ath ^ and 
rmultampk tht PAM rays and wavekngths^ 

symmetric systems. 


5 Global Coordinates 


For polarization aberrations expansions, an object and 

required. We denote the object height by H and normalize l , 

of the field of view The object is located along the global y-axis wi ou 
eene^aUty because of the rotational symmetry of the system. The entrance pupU 
roSfnaL are denoted by either (..*,) - 

coordinates In the paraxial approximation the pupils are fiat, so the z dependence 
s dropp^ eTfor convenience in much of this paper. The pupil coordinates x y and 
^ are* formalized to one at the edge of the entrance pupil. The polar angle * , 

defined so that 


x = -psin d 
y — pcos<j> . 


( 21 ) 

( 22 ) 


Figure 2 illustrates the normalized coordinate system. 

6 Paraxial Polarization Aberration Matrices 

In this section, we determine the paraxial PAM for a rotationally symmetric sur- 
face Anoendix A provides paraxial expressions for the surface normal and lo 
(T / /i btfis vectors. The single surface Jones matrix of Section 4 is simplified 
to first’ the paraxial Jones matrix and then to the paraxial PAM. The accuracy of 
the^araxial approximation for the angle of incidence and orientation of the p 

° f 1 ConsMe'r ThVpXiTain of a ray for a system of surfaces. Because the paraxial 
L-onsiaer in p (A-lO)l, we choose x-y Jones basis 

fields «»>"> " p h , h ( Jonss malM for each surface consists of a 

Shtnto t plates, the Jones matrix in s-p coordinates, and finally a 
rotation back to the x-y coordinates 


r.’i : 
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dependence of the paraxial PAM is suppressed for notational convenience The 
paraxial angle of incidence is obtained by expanding t e ang e o inci ence 

a Taylor series 


i q {H,p,<j>) = arcsin |Sj | 

1 ■ - 
6 


S„ + ; pJ + 


+2 H p cos <t> icq imq + p 2 i mj 


(28) 


.here the term lines, in object and pupil coordinate, waa retained The pa, tut- 
1,1 orientation of the angle of incidence (24) in term, of co.[2*,(£T, />,*)] and 
sin [2 0 q (H,p,4>)] is 


H 2 i 2 + 2/f pcos0»c < »m t d-p 2 4, cos2d» 
cos [2 9 q {H,p,<j>)} = • 


sin [2 6 q (H,p,<t>)] = 


-2/fpsini/i 


(29) 

(30) 


Equations (28) through (30) define the paraxial system geometry , of • ddence 

Next consider the accuracy of the paraxial orientation o P , 

and angle of incidence a, a .pherical interface. The p^aroal or.enW.on , o the 
plane of incidence uinme. that the S, and P, b„» vector, “' ID 
A measure of the accuracy of this approximation is the ratio of the intensity o 

field in the x-y plane to the total intensity 


\UA 2 + \Uy\ 2 1 

T / d '| _ J — £J ‘ = cos i 

L xy\ v ) jf/12 


(31) 


where * is the angle the ray makes with the z-axis. When l xy « 1, most of the 
light is polarized in the x-y plane and the paraxial approximation >s HJVh n 
L « 0 most of the light is polarized in either the x-z or y-z planes and the paraxi 
approximation is not valid. Now consider the angle of incidence. The surface norm 
for a sphere from the definition (A-2) is 

(32) 




where Ri is the radius of curvature of the sphere. If the axial plane waves R 
(0,0,1) are incident then 


i = arcsin 


(£) 

- i + ifZ-Y 

- R, 6 \RiJ 


1(JL X 

40 \RiJ 


+ 


(33) 


where o - ,/x 2 + y 2 Figure 3 shows a comparison between several approximations 
To he angle of incidence for a spherical mirror. The high degree of linearity of 
the exact curve, even with angles of incidence as large as 30° , permits the paraxial 
approximation to the geometry to be used for many systems. ^ ^ 
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r'ai = r a2 /r Q 

r' a4 = r a4 /r 0 - O.5(r o2 /r 0 ) 2 


(39) 

(40) 


Finally, substituting the Fresnel expansion coefficients into (25) gives the paraxial 
PAM expansion 




exp 


( a 00? + a 02q * q 

(ni2 4 


+ U04 ? + ■ • ■) <r 0+ 

,4 + . . ,)(cos 29q<r i + sin20 ? cr 2 ) 


(41) 


where », = »,(#,/>,*) and 9 q = 9 q (H,p,<P) are defined in (28) through (30) and the 
coating coefficients ao,2M and <*1,2*,? are 


do, 2 M = + 4.».« + + /2 

a l ,2t, f = [ri.jk., - r'p.jM + /2 ' 


(42) 

(43) 


The subscripts of coefficients and a ut , q are assigned as fol ^ J he fi 

subscript assumes the value 0 for the polarization 

for the linear polarization along the s-p axis. The second subscript, 2k, is the 
of the coefficient. The last subscript, q, designates surface num er. 

8 Polarization Aberration Matrices for Systems 

In this section, we obtain the paraxial PAM for an optical qfriem of ^ isotropic 
rotationally symmetric elements through fourth order. Sections 9 and 10 discus 

'^Tote^we introduce the Baker-Campbell-Hausdortf ( BCH) identity [25]. . 

Baker- Campbell- HausdorfF 1 1/ A and B are matnccs, or certain other non- 
commuting operators, then 

exp A exp B = exp [A + B + C2 + C3 + * * *] 
where C k is a linear combination of k-fold commutators of A and B, in particular 


C 2 = 2^ A ' B 1 


C 3 = l[[A,[A,B]]-[B,[A,Bjj] • 


(45) 

(46) 


The BCH identity consistently carries out operator products, retaining terms to a 

^Now^the PAM of an optical system with Q surfaces is the product of the PAMs 
of the individual surfaces 


J = Jg Jq- 1 ■ • • J 2 J 1 


(47) 


which simplifies to 

J = exp [ao<x 0 + ai<ri + a 2 <r 2 + a 3 «r 3 ] (48) 

using the paraxial surface PAM (41), the BCH identity, and the commutation rela- 
tions in Table 4. The system T? A& coefficients (to. fourth ordep are 
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J (H,p,+) 


exp < 


Poooo^o + 

Po20oH 2 <r 0 + 

Pon\Hpcos4><TQ + 

P 0022 P 2<J, 0 + 

Po 040 P 4<r O + 

Pq\ 2\H p Z COS 4> <Tq + 

Po222tfV cos 2 <t> <to + astigmatism and field curvature 
P03 ii H ^ p cos it 1 cr o 4 - distortion 

PoaqoH'vo + q uartic P' ston 

PuooH 2 <ri + 

P nn Hp{cos 4 ><Ti - sin <r 2 ) + 

PiO22P 2 (cos20 <t x - sin 2^ <r 2 ) + 

Pi 40 oP 4 ^i + 

P\z\\H z p{ 2 zos<i><r ] . -sin^<r 2 ) + 

P1222 # 2 p 2 cos 0 (cos <j> cr\ - sin <f> <r 2) + 

Pn&Hp 3 cos 4 > [(1 +cos 20 )<ri -sin 2 <£<r 2 ] + 
Pi042P 4 (cos 2 <t> <r\ - sin 2 <£ <r 2 ) + 

P3131 H P 3 s * n ^3 + circular polarization coma 
P33iiP 3 P s in0 + circular polarization distortion 
sixth order terms 


constant piston 
quadratic piston 
tilt 

defocus 

spherical aberration 
coma 


> 


( 53 ) 


or in Cartesian coordinates 


J (H } x,y) = exp < 


^ 0000^0 + 


constant piston 


P 0 200# 2 o-o + quadratic piston 

Pom^y^o + tilt 

P 0 022(x 2 + yVo + defocus 

Poo4o(x 2 + y 2 ) 2 <x 0 + spherical aberration 

Po\3\Hy{x 2 + y 2 )<r 0 +coma 

Po222# 2 (y 2 — x 2 )<r 0 + astigmatism and field curvature 
Po3iiP 3 y<xo + distortion 

P 04 ooP' 4 o’o + quartic piston 

Pl200ff 2<T l + 

P\n\H{y<rx + x<r 2 ) + 

P1022 [(y 2 - x 2 ) <xi + 2 xy<r 2 ] + 


Pl 400 -ff 4 <Xl + 

Pi 3 iitf 3 ( 2 y<Xi + X<r 2 ) + 

P 1222 /f 2 y(y<r 1 + x<T2) + 

PnaaHy [(l + y 2 - x 2 ) <n + 2 xy<r 2 : + 

Pi 042 [(y 4 “ x 4 ) <xi + 2 xy (x 2 + y 2 ) <r 2 ] + 

P 3131 Hx (x 2 + y 2 ) <X3 -feircular polarization coma 
P331 1 H 3 X&3 + circular polarization distortion 


sixth order terms 


> 


( 54 ) 


where the terms are grouped based on their H , p, <t> dependencies and P, uvw are 
the polarization aberration coefficients. The polarization aberration coefficients are 
sums over interface contributions as given in Appendix C and assigned subscripts 
using the convention: t denotes the type of polarization behavior[ 2 ], u denotes the 
order of the H dependence, v denotes the order of the p dependence, and in denotes 
the order of the <p dependence. 

The fourth order paraxial PAM paraxial in system geometry andjourth order 
in coating response to changes in'kngle of incidence. A subtle conle’quence is that 
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thin film coatings (particularly those with many layers), they can be significant. 
Our equations only calculate the system wavefront aberration contributions arising 
from coatings. The geometrical portion of the classical wavefront aberrations (those 
arising from optical path length differences) are calculated by the paraxial ray trace 

and conventional aberration calculations [18]. 

The second class of aberrations contains the amplitude or apodisaiton aberra- 
tions which are characterized by the real parts of the <ro coefficients ,JR{Pa uvwx ). 
The amplitude aberrations are variations of the amplitude of the electromagnetic 
field across the exit pupil which are independent of the incident polarization state. 
They do not describe the shape of the transmitted wavefront, only its amplitude. 
This apodisation is due to the optical system, not to intensity van a ions in the 
incident light, such as the Gaussian profile of a laser beam. The amplitude aber- 
rations describe the average of the coating amplitude transmittance of the s and 
p light (the polarization terms describe the difference). Amplitude aberrations are 
scalar aberrations and have the same functional dependence on object and pupil 
coordinates as the classical wavefront aberrations. Contours of constant apodisa- 
tion aberrations through fourth order are shown in Figure 4. Since the functional 
form is the same, the generic names of the functions have been retained with the 
prefix amplitude added: amplitude tilt, amplitude coma, amplitude spherical aber- 
ration, etc. For example, the term 3t(P 0 04o)p 4 is amplitude spherical aberration 
If mPoato) is negative the center of the pupil is brighter and the pupil becom 
dimmer quartically with pupil radius. For 3?(Po«iv«u*) positive, the pupil is ng er 
at the edge. The interpretation of all of the amplitude aberration follows m the 
same manner as amplitude spherical aberration. Intentional apodisation (versus 

apodisation aberrations) is discussed in Reference [26]. 

The third and fourth classes of aberrations contain linear diattenuation and lin- 
ear retardance aberrations and are characterized by the real and imaginary parts 
of the coefficients of <r x and <r 2 respectively. These two lasses of aberrations are 
characterized by the real and imaginary parts of the coefficients of <r x and <r 2 re- 
spectively. These two classes of aberration will be treated together under the name 
vector aberrations. Vector aberrations are conceptually different from the scalar 
aberrations since both a magnitude and orientation must be specified. The parax- 
ial vector aberration patterns through fourth order are illustrated in Figure 5. the 
length of the lines denotes the strength of the linear polarization element. The on- 
entation of the line denotes the orientation of the linear polarization element. 1 e 
patterns are the same for both linear diattenuation and linear retardance aberrations 
since both are vector aberrations. A detailed discussion of the vector aberrations is 


found in the next section. 

The effect of the vector polarization aberrations on polarized incident light is 
the same as a linear polarization element with spatially varying strength and on- 
entation Figures 6 and 7 show the effect of the three second order aberrations on 
linearly polarized light. The magnitudes of the aberrations depicted are not typi- 
cal, but have been chosen to clearly display the effect of the aberrations. Figures 
8 and 9 depict the effect of the three second order aberrations on circularly polar- 
ized light. The type and orientation of ellipse indicates the type and orientation of 
the polarization state. The position and direction of the arrow denotes phase and 
handedness of the polarization state. 

The last class of aberrations contain the circular aberrations. F igure 4 shows con- 
tours of constant circular aberrations through fourth order. The real and imaginary 
parts of the coefficient of <r 3 correspond to the circular diattenuation and circular 
retardance aberrations respectively. The circular aberrations are variations of the 
circular diattenuation and circular retardance across the exit pupil. The imaginary 
part of the term P 3131 [P3311] produces a wavefront with 3(P 3 131 ) coma P( p 33ii) 
distortion] when right circularly polarized! light, isL v pc|dent and -p'CPaiai) coma 
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Second order 




Fourth order 




Figure 5: Paraxial vector polarization aberration patterns through fourth order 


'■iti'A 


' h! 
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Linearly polarized light incident 

(a) 


•P1022 — 

(b) 


/ 


\ 


t C 

\ 


G I 
/ 


Pun — i x /2 

(c) 


Pi 200 — 

(d) 


Fieure 7' Effect of linear retardance aberrations on linearly polarized light. Par 
(a) shows uniform linearly polarized light in the entrance pupil. Parts (b) (c) 

and (d) show the polarization state across the exit pupil if the system has only the 

aberration P1022 = V r/2. Pun = J>/2, and P l20 o = J *!*■ 


if 


L * r : : 
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c % c 
o^c 

ccece 

c^e 

° e c 

O f 0 
0 ^ C ^ 0 

0 \ 0 
0 

Circularly polarized light incident 
(a) 

P1022 = J 7r /2 

(b) 

o tc-o 

- 0 A 

0 ,A 0 

0°0 , 
000 0 0 

o°fo 

0 

Pun = 

(c) 

P1200 = i ^/2 

i'i) 

Figure 9 : Effect of linear retardance aberrations on circularly polarized light. Part 
(a) shows uniform circularly polarized light in the entrance pupil. Parts (b), (c), 
and (d) show the polarization state across the exit pupil if the system has only the 
aberration Pio2‘2 = < Pun = J ir/ 2 , and P1200 — i*V 

, ' "j 
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H = 0 


H=H 


o 


H =2 H 0 


Figure 10: Paraxial angle of incidence at a spherical surface for three field an- 
gles. The length of the lines denotes the magnitude of the angle of incidence. The 
orientation of the lines denotes the orientation of the plane of incidence. 



' Figure 11: The first three principal vector aberration patterns for isotropic rota- 
tionally symmetric interfaces. The Pi 02 2 pattern varies quadratically with radius. 
The P 1042 pattern varies quartically with radius. Finally, the P 1062 pattern varies 
with radius to the sixth power. 


moves off-axis. The non-principal vector aberration patterns add to the principal 
pattern to give a decentered view of the principal patterns The aberrations change 
as the paraxial angle of incidence function changes for off-axis object points. In fact, 
the set of nth order vector aberration patterns is complete when the nth order pnn- 
cipal pattern can be decentered with a linear combination of the nth order patterns. 
Again these principal patterns are analogous to spherical aberration. In classical 
aberration theory when the object is on-axis, a rotationally symmetric interface 
only has spherical aberration (but many different orders of spherical aberration). 
As the object moves off-axis, the other aberrations are introduced. 

The completeness of the set of second order vector aberration patterns can be 
addressed by constructing a shifted second order principal pattern from a linear 
combination of the second order vector aberrations. Figure 13 shows the second 
order patterns P 102 2 , P 1200 , and P U n adding to give a decentered view of the prin- 
cipal pattern. Figure 13 (a) shows the superposition of the Pi 20 o and P l02 2 patterns. 
The vector aberration patterns (arrays of weak linear diattenuators and/or arrays 
of weak linear retarders) add as PAMs, not as vectors! Two orthogonal weak linear 
diattenuators of the same magnitude , r, add to rero diattenllatioq.a^d- reduced ampli- 
tude (transmission). Two orthogonal weak linear retarders with equal retardance 
add to zero net retardance. Figure 13 (b) shows the result of adding P 12 oo and 
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p +p 
r l200 1022 


P +P 

r 1200 1022 


/ + ~ h \ 

. -V 


p + p +/> 
r 1200 1022 1111 


p + P + P 
r 1200 r l 022 1111 


Figure 13: Addition of second order aberration patterns to give a decentered view 
of the second order principal aberration pattern 
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Table 6 : Symmetry of pupil aberration sections of rotationally symmetric systems 



On Axis Objects 

Off Axis Objects 


Tangential 

Sagittal 

Aberration type 
Aberration orientation 

Symmetry 

Linear 
Radial or 
tangential 
Rotational 

Linear 

Horizontal or 

vertical 

None 

Elliptical 

Arbitrary 

Odd about 
the y-axis 


Note: An aberration has odd symmetry if the pattern is mirror symmetric about 
some line. 


Tangential 


Sagittal 


/ 

' / 

» V 

\ V 

\ 


I ■ 


\ 

\ 

\ 1 • 

/ 

I t 

1 1 
\ '- 
\ 

• 1 

s / 

/ 

. ^ 

(a) 


« 1 • 

• 1 , 

(c) 

(d) 


\ 


/ 


\ 

I * 

/ 


I • 


Figure 15 : Example pupil aberration section for on axis object. Exit pupil aber- 
ration maps with the tangential and sagittal sections highlighted are shown in (a) 
and (b), respectively. The tangential and sagittal sections are shown in (c) and (a), 
respectively. The sagittal section is rotated 90 degrees for a more compact display. 
The aberration is P1022 - Pi 042 + P33 11 • 
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Field height Tangential section 


Sagittal section 


H=0.7 | | I I I I | 1 

H=0 | i • • • i | I 1 


Figure 17: Tangential and sagittal pupil aberrations plot for a system with nth 
order principal aberrations and piston 


can be either diattenuating or retarding giving polarization dependent apodisation 
or phase, respectively. The circular aberrations described the circular retardance 
and circular diattenuation introduced by an optical system. For systems o iso ropic 
rotationally symmetric surfaces, the circular aberrations appeared first at four 
order. Vector aberrations were conceptually different trom classical aberrations^ 
They described the linear retardance and linear diattenuation variation introduced 
by an optical system. These polarization aberrations were decomposed into a set 
of vector patterns which were each attached a weight or an aberration coefficient. 

The theory described in this paper applies to many optical systems built today. 
Thin films have a much deeper role in optical system design than merely changing 
the transmittance of a system. Thin films induce polarization aberrations or, if the 
designer is clever, control the polarization aberrations. This extension of aberration 
theory was made possible by including a Taylor series expansion of the Fresnel 
coefficients for each interface in the optical system. The resulting aberration theory 
allows the integration of thin film design and optical design for polarization critical 

optics. 
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A Paraxial Basis Vectors and Fields 

In this appendix, we derive the paraxial quantities necessary for the paraxial PAM 
in Section 6. First, general egressions for the surface n.Qftnal and local basis vectors 
are presented. Next, the paraxial approximation is applied to these quantities. The 
paraxial field is shown to lie in the x-y plane. 
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where 


Vcq _ r g sm u cq 
ftq R\q 

Vmq Jj_ 

Rq Riq 


(A-U) 

(A-12) 


and : e? and « m , are the angle of incidences for the chief and marginal ™ys respec- 
tively. This follows directly from linear nature of the paraxial raytrace [lb], ine 
paraxial system geometry is entirely contained in i c and * m . It is not necessary to 
work with the radii of the individual lenses, entrance pupils, etc. The asphencity 
and wavefront aberrations have a third order effect on the surface normal and local 
(R q ,S v P q ) basis vectors. Both 5, and P t are in the x-y plane. There is no z 

component of the electric field. 


B Polarization at Uncoated Interfaces 

This appendix examines the polarization on reflection and transmission from an 
uncoated interface. We determine the interface coefficients of Chapter 7 wit a 
Taylor series expansion of the Fresnel coefficients. The accuracy is discussed tor 

reflection by gold at A = 10.6 pm. 

The Fresnel coefficients expanded in the angle of incidence are 


1 . 2 ,V 2 - - 3 . 4 


r v = 


U = 


u = 


sin(i — i f ) 
sin (i + i 7 ) 

jvTT ( L " N l * + un* ' + 

tan (i — i') 
tan(i + i') 

N-l (. 1 , 2 . oA 2 - 6A + 9 4 , 

F+Tv 1+ iV i2iV 3 ) 

2 cos x sin(i') 
sin(i + 

2N 


N - 2 

sin (j - 
sin (i + 1 ') 

2 iV 


j\r - 1 3iV 3 + 3iV 2 - 7 N 

1 + + 24 




N + 1 


1 + 


N{N - 1), 2 , N{N - l)(9iV 2 ~ W + 5) , 4 

Z * 4" n! ' 


24 


(B-l) 


(B-2) 


(B-3) 


(B-4) 


where r , and r p are reflection coefficients, t, and t p are reflection coefficients, i is the 
angle of incidence in radians, the ' denotes quantities after the interface, N = n /» 
is the ratio of the refractive indices, and Snell s law 


n sin i = n 1 sin i' (B 5) 

was used [16,15]. The interface coefficients (B-l) and (B-2) ^exponential form are 

*■ AJ. T ^ 

- 1 ' ’ 

r,(i) = exp[r 0 + r 2 i 2 H ] (B-6) 

r p (i) = exp[r 0 - r 2 i 2 + • • •] (B-7) 
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Angle of Incidence (deg) 



s 



(b) 


Figure 18: The amplitude reflectance (a) and phase change (b) on reflection from 
a gold coating at 10 6 fim. Curves' are shown for both the s and p components ot 
the incident light. Vain s were computed using a refractive index, n = 8.8 - 64j. 
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Polarization aberrations. II. Tilted and decentered 

optical systems 


James P. M c Guire Jr. 
Russell A. Chipman 


May 21, 1990 



Many optical systems built have tilted or decentered elements. These include 
unobscured systems and systems with fold flats. Because of typically larger angles 
of incidence, polarization aberrations can be significant in these systems. Two types 
of tilted and decentered systems composed of rotationally symmetric elements are 
examined. One is systems with collinear centers of curvatures but with decentered 
pupils Symmetry in such systems allows the analysis to proceed along lines very 
similar to those in Paper I. The other is systems with arbitrary tilts and decenters. 
In these systems, the field dependences of the aberrations from each surface are not 
concentric. The extension is made by using a polarization aberration matrix with 


vector, instead of scalar, arguments. .... . , , 

The extension to tilted and decentered systems used in this paper is based on 
the principle that each surface has an axis of symmetry; and these aberrations 
can be found in the conventional fashion. Buchroeder used this principle to design 
systems composed of tilted and decentered elements [I] Thompson used vector 
algebra to combine the aberration contributions from tilted and decentered elements 
[2? Rogers extended the vector techniques to explore anamorphic and keystone 
distortion due to the tilt of the object relative to the elements in the optical system 

^ Section 1 briefly introduces the basis of aberration theory for tilted and decen- 
tered systems. Section 2 outlines the coordinate system and some ve ^or operations^ 
In Section 3 we present the PAM with vector arguments which is the i basis of the 
calculations in this paper. Section 4 examines systems with decentered pupils E ' 
ample calculations for an infrared LIDAR beam expander are given ,n Section 5. 
Section 6 explores aberrations in systems with arbitrary tilts and decenters, 
second order PAM of Section 3 is manipulated into a form convenient for summing a 
system of tilted and decentered elements. The accuracy of this method is discussed. 
A simple IR scan mirror assembly is analyzed in Section 7 References to equations 
and sections in Paper I are preceded I. 


1 Overview 

Each spherical surface and its entrance pupil form a rotationally symmetric systerm 
Thus each optical surface introduces aberrations of the same form, whether use 
in a rotationally symmetric or unsymmetric system. The . enters o t e a erra ion 
contributions of each surface are displaced due to the tilt. Figure I shows a spherical 
surface, pupil, and object. The axis of symmetry is called the local axis and connects 
the center of the pupil and the center of curvature for the surface The centra/ ay 
is defined by the center of the object and the center of the pupil for the surface^ 
Both the object field and the pupils for each surface are found by imaging the object 
and the entrance pupil through each of the surfaces prior to the surface in question. 
In rotationally symmetric systems, the local axis and the central ray coincide. In 
systems with decentered pupils, the vector aberration expansion is made about the 
line connecting the centers of curvatures of the elements. In arbitrarily tilted an 
decentered systems, the vector aberration expansion for tilted optical systems 

made about the local axis of each element. 

When the center of the pupil coincides, or nearly coincides, with the center of 
curvature of the mirror (i.e. in a Schmidt system), the local axis is either ill-defined 
or so oblique that aberration expansions about it are impractical. In these cases, 
the expansion should be made about the line connecting the object and the center 
of curvature /tte’ optical systJni IS rotationallyiymmetric about this line with the 
pupil decentered and tilted instead of the object. 
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ORIGINAL PA'rE IS 
OF POOR QUALITY 


( 3 ) 


P = ( P,<t > ) 

= p x i + p y y 

where p and 0 are the polar pupil coordinates. The pupil vector is normalized so 

that p = 1 at the edge of a circular pupil. 

Vector multiplication introduced by Thompson [2] is analogous to complex num- 

ber multiplication 

Hp = (Hp,d + <t>) W 

and gives a vector result. The square root of a vector 

n/I = (VH, Q or (VH , ^ ^ 

is a vector and follows directly from the definition of vector multiplication. The dot 
product 


H p = H x p t + HyPy 

= H p cos <f> W 

gives a scalar result. The following vector identities 

2 (.4 fl) (a 2 ^ 2 ) = A 2 (AB C 2 ) + (.-P -SC 2 ) ( 7 ) 

2 (A ■ b) (a • (?) = (A 2 • SC) + A 2 [ 3 ■ (?) ( 8 ) 

2 (.i s) (ABC 2 ) = A 2 (fl 2 C 2 ) + S J (A 2 C 2 ) (9) 

are useful in the manipulation of vector aberrations. All are easily verified by 
converting to trigonometric form. 


3 Polarization Aberration Matrices with Vector 
Arguments 

In this section, we discuss PAMs for a single surface with vector arguments. The 
PAM with scalar arguments is rewritten first for arbitrary object onen ations an 

then with vector notation. . q 

First consider the second order terms in the PAM derned in Section 1.9 


J (H,p,<i>) - exp 


Poooo^o + 

Po20oH 2 <t 0 + 

Pom^ P cos< ^ + 

P(i020P 2<r 0 + 

Pi2oaH 2 <Ti + _ • 

Pxjn,Hp{cps<j>.<ri - sin0 ( <72) + 

P\027p 7 (cos"' 2 <^ <Ti - sin 2 * to) 


constant piston 
quadratic piston 
tilt 

de focus 


(10) 


Equation (10) describes the aberrations for objects of height H located along the y- 
axis The pupil dependence is described by polar coordinates (p, *) with 4> measured 
from the object (which lies on the y-axis). If the object ,s located 0 from the y-axis, 

then the PAM is 
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Figure 3: IR LIDAR off-axis beam expander (a) and the equivalent rotationally 
symmetric system (b). In systems with decentered pupils such as a), the aber- 
rations are easily calculated by analyzing the equivalent -otationally symmetric 
system (b) and then decentering the pupil. 
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Table 1: Optical design of the LIDAR off-axis beam expander and the associated 
chief and marginal angles of incidence 


Overall 


F-number 

1 

Expansion ratio 

16:1 

Field of view 

50 prad 

Diameter 

110 cm 

Coatings 

Gold 

Wavelength 

10.6 pm 

Parabolic primary 

Focal length 

iiu cm 

Diameter 

110 cm 

Off axis section diameter 

30 cm 

Beam displacement from optical axis 

50 cm 

Parabolic secondary 

Focal length 

b.o/o cm 

Diameter 

7 . 5 cm 

Beam displacement from optical axis 

2.5 cm 

Paraxial angles of incidence 
Chief 

At primary mirror (i c i) 

—50 /irad 

At fold mirror (i C 2 ) 

50 /irad 

At secondary mirror (i c3 ) 

424 /irad 

Marginal 

At primary mirror (i m i) 

750 mrad 

At fold mirror (i m 2 ) 

-500 mrad 

At secondary mirror (ima) 

72.7 mrad 


Table 2: Aberration coefficients for the LIDAR beam expander 

Poooo = — 0.012 — 0.092.7 
P0200 = 0 
fbiu = 0 

fo020 = 9 . 

P im =3.5 x 10- 10 + 2.8 x 10" 9 j 

' ' r 'y'l'?':’ fi:!n-=''T 1 -. 3 x TO"’’ -* 1 4 x 

,T? *' ‘ p,n„ = 0.002 4- 0.013 1 
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("■') = 


exp < 


E? •Fbooo» <T o + 

constant piston 

E ? Po2oc>f (h ~ hq) ' — 

) <To + quadratic piston 

E? Pom« _ hi) ' P* 0 + 

tilt 

E, P0020lP P*0 + 

defocus 

E,*200, (H-hX A + 


Z q Pnn q ( H-h,)p A + 


E, P l022qP 2 ' A 

> 

( 


which may be rewritten 


(h,p 

^ = exp < 

P\ 200 
| Pun 1 

(ji - 01200^ “ &120oj ' A + 

7 H - a ml ) p - A + 



L Pi022P 2 * A ' 


( 15 ) 


( 16 ) 


where the polarization independent aberrations Pq U vw were dropped and the aber 
ration coefficients are 


P 1200 — ^ l200 $ 

(17) 

9 

ai200 - p l2QQ '^2 Pl200| l h,l 

q 

(18) 

£?200 = “1200 - p^Y, Puoo r'i 

q 

(19) 

Pull = y] Pullq 

q 

(20) 

ami - p nn ^ Pnnqhq 

(21) 

P 1022 = Yl Pl022q ' 

(22) 


q 


Equation (16) describes the polarization dependent aberrations through second or- 
der for systems of tilted and decentered elements. The interpretation of each aber- 
ration is discussed below. 

The first aberration to consider is polarization defocus 


P 1022 P 2 • A . ( 23 ) 

Since each contribution to the defocus terms is independent of H and therefore h, 
the aberration is independent of the tilts or decenters of the elements. 

The next term to consider is polarization tilt 


. ” aniiYp’ 'A 


(24) 


which shows the aberration resulting from a sum of tilts is of the same form as 
a single surface tilt Figure 6 shows the pupil aberration for the tilt aberration 
in (a) rotationally symmetric system and (b) tilted and decentered system. The 
aberration is centered in the field at a„„- The strength of the aberration is given 
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* 



(b) 


Figure 6 The pupil aberration map for the tilt aberration. Examples for ( a ) r °~ 
tationally symmetric systems and (b) tilted or decentered systems are shown. The 
object vectors are superimposed over the pupil to show the effect of the tilts and 
decenters on the orientation of the aberration pattern. 


11 


PRECEDING PAGE BLANK NOT FILMED 





Fieure 8' The pupil aberration mkjj'for quadratic’ piston aberration. Examples for 
(a) rotationaffy symmetric systems and (b) tilted or decentered systems are shown. 
The object vectors are superimposed over the pupil to show the effect of the tilts 
and decenters on the orientation of the aberration pattern. 
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Scan mirror 


Pupil 


(a) 




Figure 10: Infrared scan mirror system pointed at (a) 22.5°, (b) 45°, and (c) 67.0°. 
The fixed pupil of the optical system is shown. 
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Polarization Analysis of LIDARs 

Dr. James P. M c Guire, Jr. Dr. Russell A. Chipman 
Physics Department 
University of Alabama in Huntsville 
Huntsville, Alabama 35899 

June 8, 1990 


1 Introduction 


The objective of LIDAR systems is to accurately measure atmospheric : winds. The 
measurement proceeds as follows. A circularly polarized beam is sent by the L D R 
into the atmosphere. Particulates backscatter some of the light back into the LIDA 
and change the handedness of the circularly polarized light. T e . return ®‘^ 
combined with the local oscillator at the heterodyne receiver using a P^larumg 
beamsplitter. Doppler shifts and trip time of the return beam provide velocity a 

ran Any difference* in the polarization state of the return beam from the expected 
circular polarization state, decreases the fraction of the return signal combined 
with the local oscillator. Thus changes in polarization due to he 0 P U “ befo " 
beamsplitter, reduce signal and instrumental accuracy. This loss of * lgnal £ 
instrumental polarization can be minimized during the design phase, if polarizati 

18 Tlfof the mirrors in LIDAR systems change the polarization state of the light 
because the rays strike at non-normal incidence. The polarization change depen s 
on field position, object position, and incident polarization state. Polarization aber- 
rations couple some portion of the incident polarization state into orthogonal p<> 
iarization state. Polarization analysis of optical systems is reviewed by Chipman[lj. 


2 Results 

Significant contributions in two areas of polarization analysis were made, aberration 
theory and polarization raytracing. The results are described in two papers ; which 
were completed with funding from this contract[2,3] and have been submitted to 
Optical Engineering for publication. These papers are the final report for this 
contract. This contract also provided partial support for the completion i of Dr 
M c Guire’s dissertation [4], Some important aspects of this research are highlighted 

b e low i 

Polarization aberration theory describes the low order polarization in an optica^ 

svstem with a Taylor series approximation. This approximation is particularly goo 
in reflective IR systems because of the high indices of IR materials. Section II. 5 
uses polarization aberration theory to analyze a NASA IR LIDAR beam-expander. 
Polarization aberrations were found to couple less than 1.0% of the light into the 

orthogonal polarization state. . , . 

Polarization raytracing is the analogue of conventional raytracing with the po- 
larization modifying properties of the system calculated instead of optical path 
differences. Current commercial polarization raytracmg codes [5,6] are good at cal- 
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i . , i , . .. . _ svstem. However the results are not presented 

w£“o ° ZT h Sere.-nding * «» potation aberrations 

would make thadasignofUD ARotanyP > have been illili-ed , ith , en5 

^r^V5:X«r.mp,ov.d polarisation ana!,,, software should result- 
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